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^ ■ Abstract 



The low energy dynamics of a certain D-brane configuration in string theory is de- 
scribed at weak t'Hooft couphng by a non-local version of the Nambu-Jona-Lasinio model. 
^ ' We study this system at finite temperature and strong t'Hooft coupling, using the string 

theory dual. We show that for sufficiently low temperatures chiral symmetry is broken, 
^' while for temperatures larger then the critical value, it gets restored. We compute the 



latent heat and observe that the phase transition is of the first order. 
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1. Introduction 

Over the past few years there has been considerable progress in studying QCD-Uke 
theories using holography (an incomplete list of references includes |ll|-0]). The gravity 
description is useful for studying properties of the hadrons as well as qualitative features 
of QCD such as confinement and chiral symmetry breaking (xSB). The latter stands for 
the dynamical breaking of the U{Nf) x U{Nf) symmetry acting on Nf right handed and 
left handed quarks of the high energy Lagrangian down to its diagonal subgroup, U{Nf). 
In QCD the energy scales of confinement and xSB are approximately the same. 

An interesting example which exhibits chiral symmetry breaking is the Nambu-Jona- 
Lasinio (NJL) model [§], whose Lagrangian contains left and right-handed quarks interact- 
ing via a non-renormalizable quartic term. The non-renormalizability of this interaction 
makes it natural to attempt a well-defined UV completion. Recently, such a model was 
proposed within the context of string theory . The weak coupling string theory descrip- 
tion of the model considered in ||^ involves a set of Nf D8-branes and Nf anti-D8-branes 
separated by a distance L in the x'^ direction and Nc D4-branes extended in the {x^, . . . , x^) 
directions. At low energies, the system reduces to the Nf left handed quarks and Nf right 
handed quarks interacting via a five- dimensional gluon exchange. This induces a non-local 
Nambu-Jona-Lasinio-type interaction in the low energy Lagrangian for the quarks. 

An important feature of the large Nc model studied in is the separation of scales. 
Namely, the energy scale associated with confinement is well below that of chiral symmetry 
breaking. In addition, by compactifying the x"^ direction and varying the radius of the 
resulting circle one can smoothly interpolate between the NJL model and QCD. The weakly 
coupled region of the former can be analyzed within field theory |0] . At strong coupling, 
one can use holographically dual string theory. Chiral symmetry breaking and meson 
properties can be analyzed within the DBI action for the Nf D8 branes living in the 
background created by Nc D4 branes. As long as Nf « Nc, back reaction of the D8- 
brane on the background metric can be ignored. The conclusion of is that at strong 
t'Hooft coupling chiral symmetry is dynamically broken, since the corresponding brane 
configuration is energetically favored over the one with unbroken symmetry. 

In this paper we study the effects of finite temperature in the limit of strong t'Hooft 
coupling. We find that there is a phase transition at the temperature Tc ~ 0.15L~^. 
Below this temperature, there exists a thermodynamically favorable brane configuration 
where chiral symmetry is dynamically broken. As we increase the temperature above the 
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Tc, we find that the solution with the unbroken chiral symmetry is thermodynamicaUy 
preferred. As temperature reaches T^ ~ 0.17L~^ the phase with broken chiral symmetry 
ceases to exist. These properties suggest that the phase transition is of the first order. 
To confirm this we compute the latent heat, and show that it is positive. The rest of the 
paper is organized as follows. In the next section we introduce the brane model at finite 
temperature and discuss the features of the phase transition at T = Tc. We discuss our 
results in Section 3. 

Note added: As we were completing this paper we learned of the paper by O. Aharony, 
C. Sonnenschein and S. Yankielowicz [ITO] which partially overlaps with our results. 



2. D8/D4 branes at finite temperature 

We will consider the near-horizon geometry of the D4-branes in the limit studied in 
0, where the direction transverse to the D8 branes is non-compact. The near-horizon 
geometry of the L>4-branes at finite temperature is: 

ds' = (^) ' {dx^dx^ + f{U)de + {dx^f) + {^ ' (-^ + U^d^tj (2.1) 

where f{U) = 1 — Uj^/U^. Here t is the Euclidean time, z = 1, . . . , 3; t/ and O4 label the 
radial and angular directions of [x^, . . . , x^). The parameter R is given by 

R^ = T^gsNJl = nX (2.2) 

where in the last equality and in the rest of the paper we set ct' = 1. The four brane 
geometry also has a non-trivial dilaton background. 

Finite temperature implies that t is a periodic variable, 

tr^t + (3 (2.4) 



On the other hand, in order for (|2.1|) to describe a non-singular space, t must satisfy 

t = t + T7^ 2.5 



Hence, the temperature is related to the minimal value of U, denoted by Ut, as T = 

3f/^/V47ri23/2. 

Quarks can be incorporated into this system by adding Nf D8 — D8 brane pairs. 
In the following we consider Nf = 1 case for simplicity, but the discussion applies for 
finite Nf << Nc as well. At weak coupling, D8 and D8 branes are separated by the 
distance L in the a;^ direction. At strong coupling A >> L.lg, the shape of D8 — D8 
brane pair is determined by the equation of motion that follows from the DBI action. As 
we will see below, there are two possible types of solutions, straight and curved. The 
straight solution is simply r = const, and describes separated D8 and D8 branes. In this 
situation, chiral symmetry is not broken. In the Lorentzian version of (|2.1j ), the point 
U = Ut corresponds to the horizon of the black hole. The straight branes therefore cross 
the horizon in the Lorentzian space. The second type of solution is the one where the 
D8 — D8 brane surface is connected and is described by a curve U{x4) in the geometry 
( p.l|) . This configuration breaks chiral symmetry and does not cross the horizon in the 
Lorentzian space. To understand which configuration is thermodynamically preferred, it 
is necessary to compare the values of the free energy, which we do below. 

In the following it will be convenient to parameterize the D8 — D8 brane surface by 
the embedding function r(a;^, U) = x^{x^, U). The classical trajectory is x^ independent, 
and U^oo as r^ ± L/2 Q. The D-brane has two branches which are interchanged by 
r^ — r. We will mostly restrict attention to one branch, where U varies from cxd to the 
limiting value Uq at t = 0. The induced metric on the brane in the background ( p.lj ) is 

ds^ = (^] {{v^. + d^rd,T)dx^dx'' + {f{U) - l)dt^) + 2 (^^ durd^rdx^'dU 

(2.6) 
where /U = 0, . . . , 3. At this point we assume that the brane configuration is independent 
of x^. The DBI action is 

S = ^'^''^' ' j d'x j dUU''/^^l + iU/Rrf{U)iduT)^ (2.7) 

where Tg = l/(27r)^5fs is the tension of the brane and Vs = 7r^/2 is the volume of the unit 
four-sphere. It is convenient to make a change of variables 

2^3/2 ^ 4^3 



y = ^fr^ u = -Y- (2-^ 



which brings ( |2.7| ) to the form 



S = ^Ll!l^ I d'^ I dyy-\/l + fiy)idyrr (2.9) 

where 

(2.10) 

The space is restricted to he in the region y G (0, ^t), whose upper hmit corresponds to 
the black hole horizon which shrinks to a point after the Wick rotation. The equation of 
motion in the new variables is 
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l + fidyry 



(2.11) 



We denote conserved quantity in the square brackets by i/q ; this parameter parameterizes 
the solution. Then, ( |2.11|) can be written in the form 

d^J ~lWMy¥hF^ ^^'^^^ 

Note that the brane extends all the way from y = (which corresponds to the asymptotic 
region, U-^oo) to y = y^ which is determined by 

(1-4)4^-1 = (2-13) 

vt yl 

In particular, for yr >> yo^ y* ^ yo^ while in the opposite regime yr « yo, 

y. = yT(l-^ + Oiyr/yof^^ (2.14) 

In the limit yQ—>-oo the D-brane solution gets closer and closer to the horizon. 

In the discussion above yo is a constant which parameterizes the solutions of equation 
of motion. Its relation to the physical parameters yx = 3/27rT and L is determined by 

L ^ i-y* dy ^ r* dx 

2 Jo Vfiy)ifiy)yoVy''-'^) ^^ J^ ^f{x)U{x)xf/x^^-i) 

where we introduced the rescaled variable x = y/yr- Analogously, xq = yo/yr and x^ = 
y^/yr- In the limit of small yo, L differs from it by a multiplicative factor of order unity, 

.^y^mipim, ,„«,, (2.16) 



In the opposite regime we need to estimate the integral, 

L Z"^* dx 

77 ~2/T 



V(l-x6)((l-x6)a;iV2;o^-l)' 



Xq— >C)0 



:2.17) 



To estimate the behavior of the integral in (p.l7| ) at large xq it is convenient to introduce 
a new variable 



X 
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^ = (1-^ )-Ia6 



X 



omitting terms subleading in I/xq^ we obtain 
which goes rapidly to zero as yo^oo. 



Vt 



2/0 



(2.18) 
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Fig 1. The value of L/2 in the units of yr [eq. ( \2.1^ )] as a function of xq = yo/VT 



We are thus led to the following picture. The integral over x in ( p.l5| ) goes like xq for 
small Xq and peaks around xq ~ 1. It then goes down and goes rapidly to zero as xq^oo. 
Hence, for a sufficiently small temperature T « 1/L we generically have two solutions for 
a D8 brane which asymptotes to x^ = ±L/2. This picture is confirmed by the numerical 
evaluation of the integral in ( |2.15| ) shown in Fig. 1. 

For a sufficiently high temperature, there are no curved solutions, although there 
still exists a solution with constant x^. In Minkowski space, this solution describes two 
collections of branes and anti-branes which cross the horizon. For this configuration chiral 



symmetry is restored, which is what we expect from the high temperature phase. To 
understand the details of the chiral symmetry restoration, we need to compare the values 
of the action for the solutions described above. 

The difference between the values of the free energy for the straight and curved solu- 
tions is 



hr. 



straight -^curved 



T 
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straight 
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VT ry, 

dy{y-^ - 0) + / dyy-"" | 1 - ( 1 - 
y, Jo 



,16 



2/16(1 _ ^6/^6) 



-1/2N 



(2.20) 



When the temperature is small, yx » L, it is easy to consider the branch with y^ ~ L. 
There, the result is the same as in ISH, 



55 ~ 2/0 ^ ~ yj^^iyr/Lf 



:2.21) 



In particular, it is positive, i.e. the curved solution is preferred and chiral symmetry is 
unbroken. Moreover, dS blows up as L^O. 
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Fig 2. The value of 5S in appropriate units [eq. (\2.2C\ )J as a function of xq 
Vo/Vt- 



Consider now the second branch. Now the integral in ( p.20|) is dominated near the 
region y k, y^. To see this, introduce a new variable z defined by ( ^.181 ). The relevant part 
of the integral in (|2.20|) is 



(i-»/i + i: 



-16 



-Xq logxo 



(2.22) 



where we neglect terms that are 0{xq ) The result therefore is negative, and goes to zero 
as L^O: 

SS ~ -2/^"(yT/l/o)-''log(yo/yT) ~ -y^' {L/yr)^ log L (2.23) 

where we used ( |2.19| ). In fact, as we see on Fig. 2, the difference ( p.20| ) is positive for 
sufficiently small yo, but then becomes negative. 

The picture therefore is the following. When the temperature is low compared to 1/L, 
the system resides in the state where chiral symmetry is broken and yo ~ L. This state is 
described by the graph in Fig. 1 with xq « 1. From Fig. 2 it is clear that the action for 
this branch is smaller then that for the straight brane solution, where chiral symmetry is 
restored. It is interesting that there exists a second solution, with large yo for a sufficiently 
low temperature. This solution also breaks chiral symmetry, but is not thermodynamically 
preferred. 

When the temperature is raised above the critical temperature Tc ~ 0.15L~^, the 
system undergoes a phase transition: the straight brane becomes thermodynamically 
preferred, and chiral symmetry gets restored. (The critical temperature corresponds to 
Xq ~ 0.885, according to Fig. 2. This translates into L/2yT ~ 0.161, according to Fig. 
1. Using ( |2.10| ) gives the value of Tc quoted above). As the temperature is increased past 
T^ ~ 0.17L~^, the only available solution is a straight brane. This picture suggests that 
the transition is of first order. Indeed it is easy to see that the latent heat is non-zero at 
T = Tc. The simplest way to see this is to note that the entropy density jumps at the 
transition 

^'-' ^ ^\'' straight ''curved) 



V^ Fs dT 



(2.24) 



Using (|2.20|) , and computing the derivative at T = Tc we find 



The latent heat per unit volume is positive and equal to 



^^ 2.WTsVs\'(%\ T! (2.25) 



^^ ^ Tc{Sstra^ght Scurved) ^ Q.OSTT^iVcA^Tj (2.26) 



3. Discussion 



In this paper we considered the string theory description of the chiral symmetry phase 
transition. The corresponding D-brane system at weak coupling is described by a certain 



NJL model with non-local interaction of the quarks. We argued that when the system 
is strongly coupled, A >> L, the following scenario takes place. As the temperature is 
low compared to the scale set by the brane separation, 1/L, the curved D8-brane solution 
which breaks chiral symmetry is thermodynamically preferred over the straight D8-brane 
solution, with unbroken chiral symmetry. This solution goes over smoothly to the zero 
temperature solution studied in 0. In particular, the maximal value of y on the D8-brane 
trajectory, denoted y* in the paper, scales as y* ~ L for small L, i.e. the D-brane is getting 
further away from the horizon at y = |/t as L decreases. As the temperature approaches 
Tc ~ 0.15L~^ from below, the free energy of the straight brane solution approaches that of 
the curved one. At T = Tc the derivative of the free energy with respect to the temperature 
jumps, which indicates the first order phase transition. When the temperature is raised 
above T^ ~ 0.15L~^, the straight brane solution, which restores chiral symmetry becomes 
thermodynamically preferred. For temperatures higher then T^ k, 0.17L~^, the straight 
brane is the only available solution. At low temperatures there also exists a second curved 
D8-brane solution, which comes asymptotically close to the horizon as L— >0. The value of 
the action for this solution is higher then that of the other two. 

To reiterate, the scale of the temperature where chiral phase transition happens is 
set by the asymptotic separation of the D8-branes L. It is interesting that this energy 
scale is much larger then that of Hawking-Page transition which is associated with con- 
finement/deconfinement transition in field theory. To understand this better, consider the 
space ( |2.ip where the x^ direction is compactified on a circle of radius i?4 . In this case, 
the scale of Hawking-Page transition is set by yy ~ -^4- The case of i?4 ~ L corresponds 
to QCD, where both confinement and chiral symmetry breaking happen at the same scale. 
On the other hand in the case i?4 >> L, where our analysis applies, the two scales are 
widely separated, as has been also argued in 0. 

There are many open questions. It would be interesting to generalize the discussion of 
this paper for a non-zero chemical potential. There one expects non-trivial phase structure 
in the (/U,T) plane. Understanding the behavior of the mesons at finite temperature is 
another direction. 

Acknowledgements: We would like to thank E. Antonyan, M. Douglas, J. Harvey, D. 
Kutasov and A. Ryzhov for discussions and correspondence. A. P. is grateful to the Theory 
Group at the University of Chicago for hospitality. This work was supported in part by 
DOE grant DE-FIG02-96ER40949. 



References 

[1] A. Karch and E. Katz, "Adding flavor to AdS/CFT," JHEP 0206, 043 (2002) 
| arXiv:hep-th/0205236 . 



[2] D. T. Son and M. A. Stephanov, "QCD and dimensional deconstruction," Phys. Rev. 
D 69, 065020 (2004) ||arXiv:hep-ph/0304182 |. 



[3] J. Babington, J. Erdmenger, N. J. Evans, Z. Guralnik and I. Kirsch, "Chiral symmetry 
breaking and pions in non-supersymmetric gauge / gravity duals," Phys. Rev. D 69, 
066007 (2004) ||arXiv:liep-tli/ 03060 18| . 



[4] M. Kruczenski, D. Mateos, R. C. Myers and D. J. Winters, "Towards a holographic 
dual of large-N(c) QCD," JHEP 0405, 041 (2004) ||arXiv:hep-th/031127q . 



[5] T. Sakai and S. Sugimoto, "Low energy hadron physics in holographic QCD," Prog. 
Theor. Phys. 113, 843 (2005) [|arXiv:hep-th/0412141|| . 



[6] J. Erlich, E. Katz, D. T. Son and M. A. Stephanov, "QCD and a holographic model 

of hadrons," Phys. Rev. Lett. 95, 261602 (2005) ||arXiv:hep-ph/0501128|| . 
[7] T. Sakai and S. Sugimoto, "More on a holographic dual of QCD," Prog. Theor. Phys. 

114, 1083 (2006) ||arXiv:hep-th/0507073|| . 
[8] Y. Nambu and G. Jona-Lasinio, "Dynamical Model Of Elementary Particles Based 

On An Analogy With Superconductivity. I," Phys. Rev. 122, 345 (1961). 
[9] E. Antonyan, J. A. Harvey, S. Jensen and D. Kutasov, "NJL and QCD from string 

theory," |arXiv:hep-th/0604017 . 



[10] O. Aharony, J. Sonnenschein and S. Yankielowicz, "A holographic model of deconflne- 
ment and chiral symmetry restoration," arXiv:hep-th/0604161[ 



